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Starting point is a given semigroup of completely positive maps on the 2x2 matrices. This semi- 
group describes the irreversible evolution of a decaying 2-level atom. Using the integral-sum kernel 
approach jll| to quantum stochastic calculus we couple the 2-level atom to an environment, which 
in our case will be interpreted as the electromagnetic field. The irreversible time evolution of the 2- 
level atom then stems from the reversible time evolution of atom and field together. Mathematically 
0^ , speaking, we have constructed a Markov dilation || of the semigroup. 

The next step is to drive the atom by a laser and to count the photons emitted into the field by 
the decaying 2-level atom. For every possible sequence of photon counts we construct a map that 
gives the time evolution of the 2-level atom inferred by that sequence. The family of maps that 
we obtain in this way forms a so-called Davies process [Q, In his book Davies describes the 

structure of these processes, which brings us into the field of quantum trajectories Within our 
model we calculate the jump operators and we briefly describe the resulting counting process. 
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turn trajectory theory 0]. As an example we consider here the case of resonance 
fluorescence. 

Our starting point is a semigroup of transition operators {T t } t > on the algebra M 2 
of all 2 x 2-matrices. This semigroup describes the irreversible evolution of a spon- 
taneously decaying 2-level atom in the Heisenberg picture. By coupling the atom to 
a quantum noise, we construct a stationary quantum Markov process having pre- 
cisely these transition operators. If we impose the requirements that the external 
noise be a Bose field, and the quantum Markov process be minimal, then the latter 
is uniquely determined. It is called the minimal Bose dilation of (M.2,T t ,g) [PH , 
where g is the ground state of the 2-level atom. 

Since this dilation is uniquely determined, any other reversible dynamical model 
which couples (M 2 ,T t ,g) to some Bose field necessarily contains this Bose dilation 
as a subsystem. Therefore, without deriving our model from an explicit Schrodinger 
equation (by performing a Markovian limit) we may safely assume it to be a phys- 
ically correct way to describe the interaction of the 2-level atom with the electro- 
magnetic field. 

We will couple the 2-level atom to the electromagnetic field by using quantum 



I. INTRODUCTION 



In this paper we want to illustrate that quantum stochastic calculus together with 
the processes studied by Davies in his book [Q, and explained in his paper with 
Srinivas | 16|| , form a suitable mathematically rigorous framework for doing quan- 



stochastic calculus pq| , ||14j| . We use a version of quantum stochastic calculus based 
on integral-sum kernels |TTJ, [12|], ||14|| , which has the advantage that we have an 
explicit construction for the solution of the quantum stochastic differential equation 
with which we will describe the coupling of atom and field. Having this explicit con- 
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struction in our hands is important for doing the actual calculations we encounter 
later on. 

To be able to discuss resonance fluorescence we have to use a dilation where we 
have two channels in the electromagnetic field. On one of them we will put a laser 
state to drive the 2-level atom. We will call this field the forward channel and the 
other one the side channel. We will then count photons in both channels. We need 
the side channel, because there we know that all detected photons are fluorescence 
photons. In the forward channel a detected photon could just as well be coming 
directly from the laser. 

For every event that can occur in the photon counters we construct a map giving 
the evolution of the 2-level atom inferred by that event. We will see that the family 
of maps we obtain, fulfills the axioms for the processes discussed by Davies Q. We 
have constructed the Davies process of resonance fluorescence. 
Using the structure theory for Davies processes we can decompose the process 
into its trajectories ||. Within our model we calculate the expression for the jump 
operators and for the time evolution in between jumps. Note that a jump in the 
system occurs the moment we detect a photon, since our knowledge concerning the 
system changes. 

Using the above apparatus we show that the resulting counting process in the side 
channel is a so-called renewal process. 



II. THE DILATION 



Let M 2 , the algebra of 2 x 2- matrices, stand for the algebra of observables of a 2-level 
atom. On this algebra we are given a (continuous) semigroup {T t }t>o of completely 
positive maps. This semigroup describes the, generally irreversible, evolution of the 
2-level atom. Lindblad's Theorem [pl| then says that T t = exp tL where L : M 2 — > 



M2 can be written as: for A G M 2 : 

k 

L(A) = i[H, A] + J2 V*AVi - Uv*Vj, A}, (1) 
i=i 

where the Vj and H are fixed 2 x 2-matrices, H being Hermitian. In this paper 
we will restrict to the simpler case where H = and there are just two V-s. This 
means there is dissipation only into two channels, the forward channel described by 
Vf, and the side channel described by V s . We choose Vf and V s such that: 

^ = (lo)' V f = K f V > V ' = K ' V i l^/| 2 + l^| 2 = l. 

This exactly gives the time evolution for spontaneous decay to the ground state of 
the 2-level atom into two decay channels, where the decay rates are given by |ft/| 2 
and \k s \ 2 - 

We want to see this irreversible evolution of the 2-level atom as stemming from a 
reversible evolution of the atom coupled to, in this case, two decay channels in the 
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field. So let us first construct the algebra of observables for these fields. Let T be 
the symmetric Fock space over the Hilbert space L 2 (M) of square integrable wave 
functions on the real line, i.e. T := C © ©^Li L 2 (R)® sn . The electromagnetic field 
is given by creation and annihilation operators on J 7 , generating the algebra of all 
bounded operators. We need two copies of this algebra, which we denote by W/, 
the field which will be the forward channel, and W s , the field which will be the side 
channel. 

The evolution over a time t of a free field is given by the second quantization of 
the left shift, i.e. the second quantization of the operator on L 2 (K) which maps /(•) 
into /(■ + t). We denote the second quantization of this operator by St- This means 
that in the Heisenberg picture we have an evolution on W/ <S> W s mapping A into 
(S; ® S*)A(S t ® S t ) (= {S- t ® S- t )A(S t ® St)), also denoted by Ad[S t ® S t ](A). 
The presence of the atom in the fields introduces a perturbation on the evolution 
of the free fields. We let this perturbation be given by a certain family of unitary 
operators {U t }teR on C 2 <8> T ® J 7 , which will be specified later, that forms a cocycle 
with respect to the shift St® St, i.e. for all t, s G R : U t + S = (S^ s ®S- s )Ut(S s ®S s )U s . 
Given this cocycle, we let the time evolution of the atom and the fields together 
be given by the following one-parameter group {T t }teR (i.e. the evolution is now 
reversible) of *-automorphisms on M 2 <8> W/ <S> W s : for all A e M 2 <8> W/ <8> W s : 

f r4 N_ / UT^-tAStUt if t > 
H j \ S_ t U_ t AUllS t if t < ' 



The solution of the following quantum stochastic differential equation 
vides us with a cocycle of unitaries with respect to the shift: 



T5|| pro- 



dU t = {V f dA* u - V f *dA u + V s dA* Sit - V s *dA Sit - -V*Vdt}U t , U = I. (2) 

In the next section we will give an explicit construction for the solution Ut of 
this equation. It can be shown (H, ||, |]TT|, [O]) that if the cocycle satisfies 
equation (|2|) we have constructed a so-called quantum Markov dilation (M 2 <8> W/ ® 
W s , {TjtgK, id ® ® 0) of the quantum dynamical system (M 2 , {?t}t>o, 5 1 ) H> [TO , 
where is the vector state on W/ iS given by the vacuum vector. This means that the 
following dilation diagram commutes for all t > (and that the resulting quantum 
process is Markov): 



M 2 
Id®/®/ 



Tt 



(3) 



T t 



i.e. for all A e M 2 : T t (A) = (id ® ® 0) (f t (A ® J ® /)). 

Let us here look briefly in the Schrodinger picture at the above diagram. If we start 
with a state p of the 2-level atom (i.e. we are now in the upper right hand corner 
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of the diagram), then this state undergoes the following sequence of maps: 

p(-) i > p(ld ® ® </>(•)) = p ® <g> </>(•) ^ p ® <p <g> 0(f t (-)) i-> 
p ® (8) 0(f t (- <g> J <g> /)) = Tr^.^ s (T t *(p <g> (8) 0)) (•), 

i.e. p maps to Tr^.^ (T u (p ® © 0)). This means that p is first coupled to the 
two fields both in the vacuum state, then they are time evolved together and then 
there is a partial trace taken over the two fields. 

III. GUICHARDET SPACE AND INTEGRAL-SUM KERNELS 

Let us now turn to giving the explicit construction for the solution of equation (Q) . 
For this we need the Guichardet space fl [[/]] of R, which is the space of all finite 
subsets of R, i.e. Q := [J ngN f2 n , where fl n := {a C R; |cr| = n}. Let us denote by 
A n the Lebesgue measure on R n . If, for n G N, we let j n : R n — > f2 n denote the map 
that maps an n-tuple £ 2 , • ■ ■ , £ n ) into the set {ti, £ 2 , • • • , t n }, then we can define 
a measure p n on fl by: fi n (E) := ^, A n (j~ 1 (£')) for all .E in the sigma field S n of 
fi n induced by j n and the Borel sigma field of R". Now we define a measure p on 
fl such that p({0}) = 1 and p = p n on fi n . This means we have now turned the 
Guichardet space into the measure space (fl, S,p). 

The key to constructing the solution of equation (0) is to identify the symmet- 
ric Fock space T with the space of all quadratically integrable functions on the 
Guichardet space L 2 (f2,p). To see this identification note that L 2 (f2 n ,p„) is, 
in the canonical way, unitarily equivalent with the space of all quadratically in- 
tegrable functions on R" invariant under permutations of coordinates, denoted 
L|y m (R n ). It is now obvious how to identify T = C © 0^Li -^|ym(^ n ) with 
L 2 (n,p) = C©0^ =1 L 2 (fi n ,p n ). 

For every / G L 2 (R) we define the exponential vector e(f) G T in the following 
way: e(f) := 10/0 © ~^f m ® ■■■■ Note that the linear span of all ex- 

ponential vectors forms a dense subspace of T . For every / G L 2 (R) we define 
the coherent vector ip(f) to be the exponential vector of / normalised to unity, i.e. 
ip(f) = exp(— i||/|| 2 )e(/). Under the above identification of T with L 2 (fi,p), the 
exponential vector (of an / G L 2 (R)) e(f) is mapped into an element of L 2 (fl,fi) 
which we denote by n(f) and which is given by: 7r(/) : fl — > C : uj \— > Ylseu /( s )' 
where the empty product Yls&d f( s ) * s defined to be 1. We will often choose for / 
the indicator function of a certain interval / C R, which we denote by xi- This is 
the function which is 1 on / and elsewhere. 

We are now ready to start the construction of the solution U t of equation @. De- 
fine the integral-sum kernel of U t (name will become apparent in a minute) to be 
the map u t that maps four disjoint finite subsets of R, o~f, a s , Tf, r s (where / and 
s stand for "forward" and "side") to the following 2 x 2- matrix, where we write 
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<jf U <7 S U Tf U r s also as {ti, t 2 , . . . , tk} such that ti < t 2 < ■ ■ ■ < t k and k G N: 

U t ((Tf,(T a ,Tf,Ts) 

exp( 



=7r(X[o,t])(<r/ Uff s Ur/U r s ) exp(-— ^F*V)14 x 

ti 



.**-**=ly*V)V*_ a ... exp(-^V), 



where 



V f if ^ G a f 

-Vj if G r, 

V s if t,- G a s 

-V* if t ? G r s . 



Then we have the following theorem of Maassen, see [TI], [ O : 

Theorem III.l: After identifying C 2 <8> JF <8> JF with L^ 2 (^ x f2, x //), the space 
of aii square integrable functions on Q x Q with values in C 2 , the solution U t : 
L^ 2 (yi x £7, /i x ji) — > Lj. 2 (fi x Q, fj, x fj) of equation (Q) is given by: 

(U t f)(Uf,U s ) = ^ / Ut((?f,(Ts,Tf,T s )f((ujf\(Tf)UTf,(u s \(T s )UT s )dTfdT s . 



CTfGuJf ' 



Now we have an explicit expression for the time evolution TJ = Ad[i7t], where [7 t 
is given by S t U t if t > and U~{St if t < 0. The family {C/t}t e K forms a group of 
unitary operators on C 2 ® .F <8> .F describing the time evolution of the 2-level atom 
and the two fields together. Stone's Theorem says that there must be a Hamiltonian 
associated to this time evolution. This Hamiltonian has been calculated recently 

0- 

IV. THE DAVIES PROCESS 

We now return to the situation in figure [| We wish to make some changes in this 
diagram and for this we need to introduce some more notation regarding Guichardet 
spaces. Let I C M be an interval. Then the Guichardet space of I is the set 
^(-0 = [J^=o^n{f), where Q n (I) = {<j C J; \o~\ = n}. In a similar way as for 
Q, which is f2(R), we can give these sets a measure structure: (f2(f), £(/), /i). 
Given a subset E of £1(1) in the sigma field £(/), we can construct the projection 
M XE : L 2 (n if i)^L 2 (n,fx): f^XEf- 

Let I be [—t, 0), then the events in E([— t, 0)) , which we abbreviate to S t , are events 
in the output field of the atom up to time t. Remember that the evolution of a free 
field was given by the left shift and that the atom is sitting in the origin. Since 
the Guichardet space representation corresponds to the photon number picture, we 
can give concrete interpretations to the subsets in E t . For instance, the subsets 
f2 n ([— 1,0)), correspond to the events "there are n photons in the output of the 
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atom into this field up to time t" . 

Now back to the situation in figure [| Suppose we have been observing the output 
in the forward and side channel of the atom up to time t with two photon counters. 
Then we are given two events Ef and E s in £ t . Since we know the outcome of 
the measurements we have to change the time evolution of the 2-level atom, i.e. we 
have to project onto the observed events (see also This is summarized in the 
following figure: 



Mo 



Ids 



>XE f ®XE s 



E&Ef,E a 



ld®<£®0 



M 2 <g> W/ ® W s 



where we have suppressed the capital letters M in the projections. The map 
£*[E f ,E s ] : M 2 -»• M 2 : A f-> Id ® (j) ® <f>(f t (A <g> Xe } ® Xe s )) is the unnor- 
malized time evolution of the 2-level atom in the Heisenberg picture given that we 
see event Ef in the output of the forward channel and event E s in the output of the 
side channel. If we are given a state on M 2 , i.e. a 2 x 2 density matrix p, then the 
probability of seeing event Ef in the forward channel and E s in the side channel 
after t seconds of observation is given by: P* [(Ef, E s )\ = Tr(p£o[i£/, E S ](I)). 
The setting is still not complete for describing resonance fluorescence. Since we 
are not driving the atom, both the forward and the side channel fields are in the 
vacuum state, at most one photon can appear in the output. We change this by 
putting on the forward channel a coherent state with amplitude z£C, defined by: 
lz t '■ W — > C : A i — > exp(— t\z\ 2 )(^ir(zx[o,t]), Air(zx[Q,t])) ■ Note that 70 is the vacuum 
state. Putting a coherent state on the forward channel mimics a laser driving the 
atom. We have suppressed its oscillations for the sake of simplicity. Now we are 
ready to do resonance fluorescence, i.e. the diagram has changed into: 



M 9 



Id®XB / ®XB s 



M 2 <g> W f <g> W s 



£t[E f ,E a 



T t 



Mo 



Id®7* t ®7o 



M 2 ® Wf ® W s 



where the map E s ] : M 2 -> M 2 is now defined by £\[Ef, E S ](A) := Id ® 

lz t ® lo(Tt(A (gi XE f ® Xe s ))- It describes the unnormalized time evolution of the 
laser-driven atom given that we see event Ef in the output of the forward channel 
and event E s in the output of the side channel. Given a state p of the atom, the 
probability of seeing event Ef in the forward channel and E s in the side channel 
after t seconds of observation is now given by: P* [(Ef, E s )] = Tr(p£ l[Ef, E a ](I)j. 
To make the notation lighter we suppres the z in E\ in the following. 
Since L 2 (Q, S, p) <g> L 2 (Q, E, p) is canonically isomorphic to L 2 (Q x Q, S <g> S, p x p) 
we can simplify our notation even a bit further. By identifying these spaces we 
can write S t [Ef,E s ] = £ l [Ef x E s ], where the righthandside is defined by: for all 
£6SjE ( ,iGl 2) t>0: £ t [E](A) := Id ® ~? Ztfi (f t (A <g> X e)), where j Ztfi is an 



abbreviation for j Zt ® 70 . We will now study the properties of the family of maps 
we defined. 



Theorem IV. 1: The family of maps {£ t [E\\t>o,E& t ®Y, t satisfies the axioms of a 
Davies process, ^j: 

1. For all t > and E G Et ® Et, £ f [E] is completely positive. 

2. For all t > and all countable collections of disjoint sets {E n } in Et ® Et 
and for all A G M 2 : £* 



3. For all t > we have £ l 



n([-t,o)) xd(ko))] (/) = /. 
fi([-t,o)) x fi([-t,o))l(A) = a. 



4. For all A G M 2 : hm t _ £* 

5. For all t,s>0 and E G E s <g> E s , F G E t <g> E t and all A e M 2 we have: 
£*[F] o £ S [F](A) = £ s+t [F - sUF](A), 

where F — s G E([— £ — s, — s) Cg> E([— £ — s, — s)) is given by: 
F - s = {(// -s,f s - s); (//, / s ) G F)} and is dehned by: 
AJF = {(cj/ U 07, u; s U cr s ); (a;/, G A, (07, cr s ) G F}. 

Proof. The only point where there is really something to prove is point 5. Let us 
first introduce some short notation which we shall only use in this proof. Let 7r(z t , 0) 
denote n(zx[o,t]) ® ^(O) and denote St <8> St just by St- Further we use the notation 
a t {U s ) for sJ t U s S t . Then for all A G M 2 , s, t > 0, F G E s ® E s and F G E* ® E 4 we 
have: 



£'[F] ( (tt(z S) 0), T S (A <g> x e )tt(z s , 0)) ) exp(£|^| 2 ) 



£*[F] o £ S [E}(A) 
exp(-(s + £)|z| 2 ) 

7r(zt, 0), f t ((7r(z a , 0), f s (A (8) x^M*., 0)) ® Xf)<z u 0) / 
7r(«t, 0), U;(n(z s , 0),f 8 (A ® x^M^, 0)) ® XF+tU t ir(zt, 0) ' 
7r(«t, 0), U;{S-tn{z s , 0),S- t f.(A ® XE)StS- t Tc(z„ 0)> ® x^E/^, 0)^ 
7r(«t, 0), F/;(S_ t 7r(z s , 0),(j t ([/ s )M ® x JW+5 <r t (£/ a ),S r _ t 7r(;* SJ 0)) ® Xf+*^(^, 0) 



Now we use the cocycle identity and the continuous tensor product structure of the 



symmetric Fock space to obtain: 
£*[F) o£ s [E]{A) 



(_ f A, , 2i ^(z t+s ,0),(at(U s )Ut)*A^>XF+tOE+t+s^t(U s )U t 7r(zt + s,0) 

6Xp(^ ~r bj\Z\ 

7r(z t+s , 0), f/; +s A <g> XF+tus+t+5^+ s 7r(^+ s , 0) 

£'+t[p - sOE](A) 



n(z t+s , 0), T t+S (A g> Xf- s us)t!"(^+ s , 0) 



exp(— (s + £)|z| 2 ) 



□ 

Define maps 1^ : M 2 — > M 2 : A i— > £* [{(0, 0)}1 (A). They represent the evolution 
of the atom when it is observed that no photons entered the decay channels. Then 
we have that for alU, s > : Y t Y s = £ t [{(0, 0)}] o £ s [{(0, 0)}] = £ t+s [{(0, 0)} - 
sO{(0,0)}] = £' +s [{(0,0)}] = Y t+S , i.e. the family {Y t } t > forms a semigroup. 
Now observe that for A e M 2 and t > we have: 

y t (A) = [({0}, {0})] (A) = id ® 7 * ® io(f t (A ® Xm ® Xm )) = 

(n(z t ) ® tt(0), ® X{0} ® X{0 } C/ t 7r(^) <8> tt(0)) exp(-t|,2| 2 ) = 
(t^) ® tt(O))*(0, 0)A(?7 t 7r(z t ) ® tt(0)) (0, 0) ex V (-t\z\ 2 ). 

If we define B t : C 2 -> C 2 : u ^ (exp(-§t|2| 2 ) (^yr(z t ) ® tt(0))) (0, 0)v, then we 
see, using Theorem |111.1| , that B t is the following semigroup of contractions: 

B t = exp ( - i(H 2 / 2 + y*F + 2zV;)t), (4) 

and for all A e M 2 : It (A) = B^ABt. We say that the Davies process £* is ideal, 
see [@. 

Furthermore Theorem [IV. 1| point 2, leads to: 

0) x O)\{(0, 0)}] (/) = /- 

If we use this and the expression for B t then we can do some estimations which 
in the end lead to: 

£*[fi[-t,0) x fi[-t,O)\{(0,0)}](7) < tKI, 

with X = (2|z| 2 |kj| 2 + 1). This property can be summarized by saying that the 
Davies process £* has bounded interaction rate, see M. 

V. QUANTUM TRAJECTORIES 

In the seventies Davies studied the structure of what we now call Davies processes 
[|J. Let us first state his results, as far as relevant, in the context of the process we 
are studying. 
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Lemma V.l: (Davies ^J) Given any ideal Davies process £ l with bounded in- 
teraction rate, as defined in the previous section, we have existence of the following 
limits: 

t 7 / :=lim-£ t [fi 1 [-t,O),{0}] and J s := lim ~£*[{0}, fii[-t, 0)] . 

These completely positive maps represent the action we have to apply on the 2- 
level atom the moment we see one photon appear in the forward, respectively side 
channel. They are the jump operations for these channels. We will explicitly calculate 
these limits later on, but first we turn our attention to decomposing the Davies 
process into its trajectories 0. For this we need the following definition. 

Definition V.2: Let Y t : M 2 — > M 2 be the maps from the previous section, i.e. 
Y t = £*[{0},{0}] and let Jf and J s be the maps from lemma |V.1| . Let ojf and 
uj s be disjoint elements of Q[— t, 0) and denote ojf U u s also as {ti, . . . ,tk} where 
— t < t\ < t2 < ■ ■ ■ < tk < for a k G N. Then we define: 

W Y , Jf ,jM^s) ■= Y tl+t J ll Y t2 . tl J^ . ..^Y_ tk , 

where J tl denotes J s if ti G uo s and J'f if tj G ojf. 

Since Y t is the time evolution of the system when, both in the forward and the 
side channels, no photons are detected and Jf and J s are the jump operations that 
we have to apply when a photon in the corresponding channels appears, it is clear 
that the string of maps Y t ^ tx J ix Y tl ^J t2 ■ ■ ■ J tk Y tk represents the trajectory of an 
observable A in M 2 when we find the outcomes ojf in the forward and u s in the 
side channel during our counting experiment. The following theorem of Davies || 
shows how to decompose the Davies process into its trajectories. 

Theorem V.3: (Davies jj^) Given any ideal Davies process £ l with bounded 
interaction rate, as defined in the previous section, we have for allt > 0, Ef, E s G S t 
and Ae M 2 : 

£ t [E fl E s ](A)= [ Wy Jf , Js (u;f,u Js )(A)d f i(u J f)d t i(u s ). 

J EfXE s 

In the previous section we already found the expression for the time evolution in 
between jumps: Y t . We now turn to the calculation of Jf and J s . For all A in M 2 
we have: 

JM) = I- ^' [PM-t, o). {9})] (A) = , im $ Ad 8 1(0)(W ' « {A)i ° 



no t 1 J tj.0 fexp(— t\z\ 
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Now look at U t Tr(z) <g> 7r(0)({s}, 0), we find by using Theorem 
^7r(z)®7r(O)({s},0)= V / 0,^0)^-^+1^ = ^(0,0,0,0) + 



aC{s} 



z 2 / u t (0,0,{r},0)dr + u t ({s}, 0, 0, 0) + z / u t ({s},0, {r},0)dr + 



o 







z 



2 



Therefore we get, for all ^4 e M 2 : 

7M , r o f Ad[^)^vr(O)({ g } ; 0)](^ 

" S texp(-t|z|2) " M 



2 



(4) 



Ad[z/ 2 + Vf](A). 
Let us now turn to the calculation of J s . We find for all A 6 M 2 : 



Taking a closer look at U t Tr(z) ® tt(O)(0, {s}), applying Theorem [LII.1 : 
[/ 4 7r(z)®7r(O)(0,{ S })= / « t (0,{ S },r,0)^ r ldr = ^(0,{ S },0,0) + 



f M i (0,{s},{r},0)rfr + 2 2 / / u t ($,{s},{n,r 2 },$)dr 1 dr 




1 ^exp(-f) 



Therefore we get, for all A £ M 2 : 

/>d^ 7 r(,)^ 7 r(O)(0,{ g })](AK, _ 
Js(A) " tff texp(-^P) " Ad 

= Ad[Vi](A). 





K. 



Since we are driving the atom with a laser now, the time evolution when we 
do not observe the side channel nor the forward channel is now given by T 4 Z := 
S f [n[-t, 0),fi[-t, 0)] and no longer by T t . We will now derive the Master equation 
for this new time evolution. For this we need the Dyson series: let Lq and J be 
maps from M 2 — > M 2 , then for all t > 0: 



exp (t(L + J)) = / exp ((c^i + t)L ) Jexp ((u; 2 — ^i)L ) J . . . Jexp(— uj k L )dui, 
Jn[-t,o) 
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where we have written u) as {a>i, . . . , u>k} with — t < io\ < . . . < uo^ < 0. 

Now remember that {Y t } t > is a semigroup so we can write Y t = exp(tL ). Then, 

using Theorem |V.3| and twice the Dyson series, we see that: 



T* = e t [to[-t,0)M-t,0)]= [ W Y ,j s ^f^s)duj f du Js 

J n[-t,o)xn[-t,o) 

exp (t(L + J f + J 8 )). 



This means we get the following Master equation: 

^-T t z = L + Jf + J s = ~{V*V,-} + [zV*-zV, -}+V* ■ V, (5) 

which is exactly the Master equation for resonance fluorescence (see ||) if we take 
z = —i^ with Q, the Rabi frequency, real. 

In the quantum optics literature (see for instance @), usually there is no photon 
counting measurement done in the forward channel, i.e. Ef = Q[— t, 0). From 
here on we will do the same, we define for all t > and E s G E 4 : .M*[i? s ] : = 
£ t t, 0), E s ] . In the following we will also suppress the index s on E s . Using 



the Dyson series and Theorem V.3 we find, for all t > and £gS ( : 



M'[E] = / W z>Js (ou)Mu), (6) 



E 

where the time evolution in between side-channel-jumps Z t is given by 
Z t = exp (t(L + J7/)) and Wz,j s is defined in the obvious way analogous to 
Definition |V.2| . Note that we have found exactly the same jump operator and 
time evolution in between jumps as in the usual quantum optics literature, see 
for instance 0, 0, i.e. we have succeeded in constructing the Davies process of 
resonance fluorescence with quantum stochastic calculus. 



VI. A RENEWAL PROCESS 

We will now look briefly at some features of the process M 1 we obtained. It is 
easily seen from the fact that {J s ) 2 = (i.e. 172(0) = 0) that the photons in the side 
channel arrive anti-bunched: the probability to see two photons immediately after 
each other is 0. We will now show that the photon counting process in the side 
channel is a so-called renewal process. 

We denote E* := E[0, t) and, via a shift, we let events E in E* correspond to events 
E — t in the output sigma field E 4 . This means that an element u> = {u\, . . . , Uk} in 
E G E* with < u\ < . . . < uJk < t should be interpreted as seeing the first photon 
appear in the side channel at time u\, the second at time 0J2 up to the A;'th photon 
at time u^. 

Given that we start the photon counting measurement in the initial state p, we 
define on the sigma fields E* (t > 0) probability measures in the natural way: for 
E G E* : P* [E] := Tr(p.M'[£ - £](/)). The family of sigma fields {E*} t > generates 
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a sigma-field S°° of Q[0, oo). Using that T*(I) = I, see equation (|), we find for all 
E G X 4 : 

P* +s [£] = Tr(p_M i+s [(£Uft[M + s)) -(t + s )](/)) = 
Tr(pM t+s [E - (t + s)Ofl[-s,0)](I)) = Tr (p.M'[£ - t]M s [fl[-s, 0)] (/) 



T^pA**^ - t]T/(/)) = Tr(p.M<[£ - *](/)) = Pp[-B], 

so P p +S [.E] does not depend on s. Therefore the family {P*} t >o on the sigma-fields 
{£*} 4 >o is consistent, hence by Kolmogorov's extension theorem it extends to a 
single probability measure P p on E°°. 

In the following, when we write u <G f2[0, oo) as {o>i,a> 2 , • • • }, we imply that < 
u)\ < LO2 < For j — 1, 2, . . . we define random variables: 

X,: fi[0,oo)^M : w = {wi,w 2 , ...}■-► | ^ J ^ 

where we take uj to be 0. These random variables give the time elapsed between 
the (j — l)th and jth detection of a photon. To prove that the counting pro- 
cess is a (modified) renewal process we have to show that for i = 1,2,... the 
random variables Xi are independent and for i = 2,3,... they are identically 
distributed. This means we have to show that for % = 2,3,... the distribu- 
tion functions F Xi (x) '■= Pp[-^i < x] are equal and for i,j = 1,2,... the joint 
probability distribution function F XiyX ,(x,y) '■= P P [^ < x A Xj < y] factorizes: 
F Xi ,Xj {x, y) = F Xl (x)F Xj (y). 

Let us first introduce some convenient notation. Note that, using equation (^), we 
have for all i? e S*: 



F P [E] = P*[£] = Tr(p^ W z ,j 3 (uj)d»(u)(I) y 
Tr(p / Z^JgZ^-^Js . . . J s Z t - Uk (I)dfJ,(u)j. 

J E 

We will denote: X\ :— x 2 := u 2 — uji, ■ ■ ■ , Xk+i :—t — u k , then we can write: 
¥ p [E) = J Tr (j)Z xl J s Z X2 J s . . . J s Z Xk+1 (I)y^u). 

Let P denote the matrix ( n ^ V then we have: 



^ fe+1 (/)= ( |Ks|2( V (/))22 0) -\^\\z Xk+1 {i)) 2 ^ 

JsZ Xk {P) = ^ 2 ( Z f P))22 0) = \K S \ 2 (Z Xk (P)) 22 P,... ., 

J s z X2 (P) = (\ K '\ 2 ( Z *W)* °) = \ K f(Z X2 (P)) 22 P 
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Therefore, if we define z(x) : = \k s \ 2 (Z x (P)) 22 , zi ast (x) := \k s \ 2 (Z x (I)) 22 and 
Zfirst(x) '■— Tr(pZ x (P)), we can write (see also jj): 

F p [E] = Zf irat (xi)^ Y\_ z ( x i)j z iast(x k +i)dn(uj). (7) 

^ E 1=2 

We would like to stress that this formula is only valid for events E G E* and not 
for all events in E°°. 

For t > we introduce the following random variables: 

N t : fi[0,oo)^N: u ^ |w n [0, t}\, 

counting the number of photons appearing in the side channel up to time t. Since, for 
strictly positive driving field strengths, i.e. \z\ > 0, the eigenvalues of the generator 
L$+Jf of the semigroup Z t all have strictly negative real parts, we have lim^oo Z t = 
0. Using this, formula (0) and the fact that the event [Nt = 0] is an element of E*, 
we obtain: 

lim F p [N t = 0] = lim z first {t) = 0. 

Now suppose we have that lim^oo P p [iV t < n] = for a certain n G N. For s < t 
we use: F p [N t < n + 1] = F p [N t <n + 1\N S < n}F p [N s < n] + F p [N t < n + 1\N S > 
n]F p [N s > n]. Therefore we have: 

lim F p [N t < n + 1] = lim lim F p [N t < n + 1] = 

lim lim (w p [N t <n + l\N s < n]F p [N s < n] + F p [N t < n + 1\N S > n}F p [N s > n]) = 
lim lim P*[JV t < n + 1\N S > n] = lim lim zi ast {t - s) = 0. 



Now using induction, we get for n G N: 

lim F p [N t < n] = 0. 

We are now ready to calculate the distribution functions Fjq and Fx i) x y The 
problem is that for instance the event \X% < x] G E°° is not an element of E* for a 
t G R. We solve this by conditioning on the event [Nt > i] and taking the limit for 
t to infinity: 

F Xi {x) =F P [X { <x] = 

lim (p p [Xi < x\N t > i]F p [N t >i}+ F p [Xi < x\N t < i]F p [N t < i}) = 

t — >OQ \ / 

lim P*K <xAiV t >z]. 
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Now we use again formula (|7|) to obtain for i > 2: 




oo „ k 



Then we use dominated convergence to interchange the limit and the integral to 
obtain: 



When i = 1 we can repeat the whole calculation to find the same result when for 
z we substitute Zfi TSt . It is now obvious that for % = 2, 3, . . . the random variables 
X{ are identically distributed. 

In a similar fashion, only extracting two integrals now, we find that for 
i,j = 2,3,... : F Xl ,x 3 (x,y) = J* J" z(x')z(y')dx'dy' . If % or j is 1 we again 
have to substitute Zf irst for z. It is now obvious that the random variables Xi and 
Xj are independent. We conclude that the family of random variables {^Q}i=i,2,... 
is a (modified) renewal process. 
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